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INTRODUCTION 

The method of sum rules is 35 years old. In spite of this respectable age, the method is being permanently enriched by 
new ideas and new calculations and remains one of the widely used and competitive tools both for the determinations 
of the fundamental QCD parameters (e.g., quark masses and a s ) and for the calculation of hadron properties. In this 
talk we review the recent progress in the applications of QCD sum rules to hadron properties with the emphasis on 
the selected topics: (i) sum rules for two-point vacuum correlation functions and leptonic decay constants of heavy 
mesons; (ii) sum rules for three-point vacuum correlation functions, form factors and three-meson couplings; (iii) sum 
rules for exotic tetraquark states. 

QCD sum rules [1] (see also [2, 3] for further references) is one of the main analytic methods for the study of 
hadron properties from the field-theoretic Green functions (correlators) in full QCD. The correlators are calculated by 
means of the Wilsonian operator product expansion (OPE) which provides the rigorous framework for the separation 
of long and short distances, in QCD being dominated by nonperturbative and perturbative physics, respectively [4]. 
The OPE clearly identifies, e.g., the origin of chiral symmetry breaking and the emergence of hadron masses, leads to 
factorization of complicated amplitudes of hadron interactions at large momentum transfers. 

• QCD sum rules provide hadron amplitudes which satisfy all rigorous properties imposed by perturbative QCD and, 
at the same time, contain nonperturbative contributions determined in a unique way. As an OPE-based method, QCD 
sum rules are formulated in the Euclidean region. However, by combining OPE with the knowledge of the analytic 
structure of the Green functions and resummation schemes, the analytic continuation to the Minkowski space may be 
performed. In this respect QCD sum rules may have a broader range of applicability than lattice QCD. Last but not 
least, as an analytic method, QCD sum rules provide physics insights in the hadron structure, which are not easy to get 
from the numerical results of lattice QCD. 

• The method of QCD sum rules favourably compares with other analytic methods, such as effective theories or 
functional methods: the method of sum rules is based on the Wilsonian OPE in full QCD and therefore involves no 
other implicit assumptions often present in other analytic method. 

a. OPE and the sum rule for the correlator 

The basic object in the method of QCD sum rules - as well as in lattice QCD - is the vacuum-to-vacuum correlator, i.e., 
the vacuum average of the T -product of quark and gluon currents. In lattice QCD, one finds this correlator numerically 
at large values of the Euclidean time T. In the method of QCD sum rules, one calculates the correlator analytically as 
the Taylor expansion in T. Technically, one considers a so-called Borelized correlator, i.e. applies the Borel transform 
to the Feynman diagrams, written as spectral representations in the energy variables. The inverse Borel mass parameter 
is related to T. The OPE provides the analytic double expansion of this correlator in form of a perturbatively calculable 
power series in the strong coupling constant a s and in powers of t; the “power corrections” — terms involving powers 
of t — are given via condensates, expectation values of gauge-invariant operators over the physical vacuum in QCD; 
these condensates describe in an unambiguous way nonperturbative QCD contributions. 

Alternatively, one may derive a representation for the Borelized correlator in terms of the intermediate hadron states. 


The two representations for the Borelized correlator — by OPE and by sum over hadron states — constitute the two 
sides of the QCD sum rule. 

b. Isolating the ground-state contribution from the Borelized correlator 

At large T, the ground-state dominates the correlator which thus fully determines the ground-state parameters. In the 
region of small and intermediate t, where the truncated OPE gives a good description of the correlator, excited hadronic 
states give sizeable contributions. In order to get rid of the excited states and to isolate the ground-state contribution 
from the correlator, one invokes the idea of quark-hadron duality [5-7]: the excited states are dual to high-energy 
parts of Feynman diagrams of perturbative QCD. The ground-state contribution is then equal to the “dual correlator” 
- the correlator in which the spectral integrals for perturbation theory diagrams are cut at a certain effective continuum 
threshold s e ff- or simply “effective threshold". The effective continuum threshold differs from the physical continuum 
threshold determined by masses of low-lying hadrons. Obviously, apart from a truncated OPE for a correlator, the 
effective continuum threshold is a crucial ingredient of every sum-rule extraction of ground-state parameters; this 
quantity governs the accuracy of the quark-hadron duality and determines to large extent the numerical value of the 
extracted parameters of the bound state. The truncated OPE itself cannot provide precise values of the ground-state 
parameters. Therefore, the method of QCD sum rules provides hadron parameters with some uncertainty which is 
referred to as systematic uncertainty [ 8 ], 

Understanding the properties of the effective continuum threshold and finding a criterion for fixing this quantity is 
the key to obtaining reliable hadron parameters from sum rules. 

1. TWO-POINT CORRELATION FUNCTION AND THE OPE 

Let us start with the simplest object - the two-point correlation function; the perturbative expansion for this object 
is known to a higher accuracy compared to more complicated correlators. Because of that, the formulation and 
application of the appropriate and reliable algorithms for the extraction of the hadron parameters from this correlator 
is becoming increasingly important. 

The two-point function, i.e. the vacuum average of the T -product of two interpolating quark currents is the basic 
object for the sum-rule calculation of the decay constants of the heavy-light mesons such as B , B s , D, D s or their vector 
analogues. For instance, for heavy-light pseudoscalar currents 75 = m^qiyff) (here TTif, is the scale-dependent MS mass 
of the heavy quark and Mb will denote its pole mass; the light-quark mass is neglected) one obtains 

n(p 2 ) = i J d 4 xe ipx (a | T ( 75 (x)jl ( 0 ))|ft) ( 1 ) 

The Wilson OPE for the T -product and for the correlation function has the following form: 

T ( 75 (x)j\ (0)) =Co(x\p)l +Y J Cn(x 2 1 p ): 0{x = 0,/i): (2) 
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Here the physical QCD vacuum |£2) is a complicated object which differs from perturbative QCD vacuum |0). The 
properties of the physical vacuum are characterized by the condensates - the nonzero expectation values of gauge- 
invariant operators over this physical vacuum: 

<fl|:d(0,/x):|fl)^0. (4) 

The numerical estimates for the condensates may be found in [2, 3]. Here we only list the recent determinations of the 
lowest-dimension condensates which claim an extremely high accuracy: 

(Cl\qq(2 GeV|£2)^ = (282±2 MeV) 3 [9], (Cl\—G a uv G a ’ pv \Q.) = 0.013 ±0.0016GeV 4 [10], (5) 

n F 

The two-point function satisfies the dispersion representation (which requires subtractions not shown here) 

U (P 2 ) = f 


(6) 



and may be calculated both using OPE (which gives it in the form E1 ope(p 2 )) and using the sum over the hadron 
intermediate states (which gives it in the form Ilhadr(p 2 ))- The sum rule is the statement that both forms represent the 
same quantity and thus should be equal to each other 

TIope {p 2 ) = Tlhadr {p 2 ) • (7) 


The spectral densities for the two representations read 


Pope(s) 


p pert (^Ai) + ^C„5W( S -M 2 )(0|(9 ; ,(Ai)|0) , 

n 
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Here Mg denotes the heavy-meson mass, fg is its decay constant defined as 

(0\j 5 \B)=f B M 2 B . 


( 8 ) 

(9) 


The truncated OPE series has quark and gluon singularities and does not have the hadron ones; therefore, comparison 
of the truncated OPE and the hadron representation in (7) may be done in the region of p 1 far from hadron thresholds 
and resonances. 

Performing the Borel transform which serves several purposes (suppressing the contribution of the excited states, 
killing the subtraction terms in the dispersion representation for n(p 2 ), improving the convergence of the perturbative 
expansion [1]) one arrives at the Borel image of the two-point function 


n h adr(T) = J dsexp{-sx)p hadr (s) = f§M B e M S X + J dse n p hadr (s), (10) 

•Sphys 

where s p h ys = [Mg* +Mp) 2 is the physical continuum threshold, determined by the masses of hadrons which may 
appear as the intermediate states, and 


rioPEO) — J dsQX p( ST)Pope(s) — J dse P P ert(SifO T IT p ower(TP) ; 


where power corrections n p0W er(T, jit) are given via the condensates and radiative corrections to them. 
The sum rule now takes the form 


( 11 ) 


EIoPeO) = Tlhadr (r)- 


( 12 ) 


Recall that the hadron (i.e. full-QCD) representation Tl| la di (T) is an infinite sum of the exponential terms, whereas 
power corrections in IIope(t) contain polynomials in T multiplied by exp(— M%x). Therefore the truncated OPE 
provides a good description of n ha dr(T) at “not too large” values of t. This determines the choice of the Borel window- 
the working T-range where the OPE gives an accurate description of the exact correlator (i.e., all higher-order radiative 
and power corrections are under control) and at the same time the ground state gives a “sizable” contribution to the 
correlator. 

The best-known 3-loop calculations of the perturbative spectral density [11] have been performed in form of an 
expansion in terms of the MS strong coupling « s (p) and the pole mass Mg: 


Ppen(s,H) =p {0 \s,Mb) + 


a s (p) 

K 


P W (sM) + 



(13) 


An alternative option [12] is to reorganize the perturbative expansion in terms of the running MS mass, mg(v), by 
substituting Mg in the spectral densities p( l \s,M%) via its perturbative expansion in terms of the running mass mg(v) 


M b 


m b (v) 


l + 



(14) 


As noticed in [12, 13], two different scales, p and v, naturally emerge when reorganizing the perturbative expansion 
from the pole ft-quark mass to the running /?-quark mass. In our discussion we do not distinguish between these scales, 
but in practical calculations the scales have been treated independently. 



Advanced algorithms for an isolation of the ground-state contribution 


The hadron representation contains the sum over all hadron intermediate states, whereas we are primarily interested 
in the ground state contribution. To exclude the excited-state contributions, one adopts the duality Ansatz : all con¬ 
tributions of excited states are counterbalanced by the perturbative contribution above an effective continuum thresh¬ 
old, s e s(x,/J.) which differs from the physical continuum threshold. Applying the duality assumption yields: 

Seff(T,M) 

ffM B e B = J dse Ppert(s,/r)+IIp 0 wer (T,/T) = II dua i(T,s e ff(T,/T)). ( 15 ) 
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Therhs is the dual correlatorYl^ ua \{x,s e s(x)) (we shall not explicitly write p as an argument of s e g but this dependence 
should be kept in mind). Obviously, even if the QCD inputs Ppert(s,jtl) and n p o W er('t : ,p) are known, the extraction of 
the decay constant requires s e ff(r,ju). Let us emphasize, that the effective threshold should be the function of x and 
p: (i) one can easily check that s e ff should depend on x in order the T-dependences of the r.h.s. and the l.h.s. of (15) 
match each other; (ii) since the truncated OPE is used in the r.h.s. of (15), the effective threshold also depends on the 
choice of the scale p. 

In early applications of the method of sum rules, it was common to use the approximation s e ff(x) = const; the value 
of this constant has been fixed by requiring the maximal stability (i.e. the least unphysical dependence of the hadron 
observable on the Borel parameter x). This procedure proved to work reasonably well, although it did not allow one 
to probe the uncertainty of the extracted hadron parameter induced by using the approximation of a constant effective 
continuum threshold. 

It should be emphasized that even if the OPE for the correlation function is known with very high accuracy in the 
Borel window, the hadron parameters can still be determined with some uncertainty which reflects the limited intrinsic 
accuracy of the method of sum rules. We refer to the corresponding uncertainty as to the systematic uncertainty. The 
latter is related to the adopted prescription for fixing the effective continuum threshold s e ff(T). 

As the accuracy of the OPE for the correlation functions has increased, one faced the acute necessity to provide 
more accurate and reliable procedures for the extraction of hadron parameters: gaining control over the systematic 
uncertainties has become mandatory [8]. 

The results of [14] demonstrated that in those cases where the bound-state mass M B is known, one can use it and 
improve the accuracy of the decay constant. We introduce the dual invariant mass M dua i and the dual decay constant 

/dual 


M dualO) = -^l0gn dua i(T,S e ff(T)), fduaM = M B 4 ^ n dualO, S e ff (?))• (16) 

The deviation of M dua i(t) from M B measures the contamination of the dual correlator by excited states. 

Starting with any trial function for .v e ir (T) and minimizing the deviation of M dua i from M B in the T-window yields a 
variational solution for .? e ff (t). As soon as the latter is found, one readily obtains the corresponding decay constant f B 
from (15). 

We consider polynomials in x and obtain their parameters by minimizing the squared difference between A7j u;i and 
M b in the T-window: 


X 2 = E [ M dual(L') -Ml] 2 . ( 17 ) 

1=1 

As shown in several exactly solvable models, the band of the estiamtes for f B corresponding to the variational solutions 
for linear, quadratic, and cubic trial s e ff(r), provides a realistic estimate for the systematic uncertainty of the decay 
constant [ 15 , 16 ]. 

The resulting f B obtained according to the procedure described above is sensitive to the input values of all the 
OPE parameters (quark masses, a s , the condensates) which are known with some uncertainties thus yielding the 
OPE-related uncertainty of f B . To obtain the latter, one assumes the Gaussian distributions for the OPE parameters 
mentioned above. Moreover, because of the truncation of the OPE series, the decay constants exhibit an unphysical 
dependence on the precise value of the renormalization scales /I. A priori, any choice of the scale is equivalently good; 
therefore, we average over the scale in some intervals assuming the uniform distribution of t/. 



• Another simple algorithm for fixing the T-dependent effective threshold in the Borel sum rule has been recently 
adopted in [17]: for each value of T the authors calculated /Vf ( | ua | (r) neglecting the T-dependence of s e ff(T) and then 
easily obtain ,s e rr by solving the equation A7 dlla | it) = Mb- Obviously, the resulting effective thresholds do depend on t; 
neglecting their T-dependence while calculating the dual mass leads to some intrinsic inconsistencies. Following our 
old idea, we tested the algorithm of [17] in a quantum-mechanical potential model for the case of a potential which 
contains the confining and the Coulomb parts [15]. This analysis shows that in quantum mechanics the algorithm with 
the variational solutions described above provides more reliable and accurate estimates for the decay constants of the 
heavy-light mesons compared with the algorithm of [17]. 

• An interesting approach to the extraction of the ground-state parameters within the finite-energy sum rule has been 
formulated and applied to the decay constants of heavy-light and heavy-heavy mesons in [18]. We have also tested 
this algorithm in the potential model [15]. For the potential-model parameters appropriate for for heavy-light mesons 
the algorithm of [18] was shown to provide rather accurate estimates for the decay constants such that the “invisible” 
systematic error remains at a few percent level only. 


Charm sector 

For the extraction of the decay constants of the charmed pseudoscalar and vector mesons, one makes use of the 
best-known three-loop expression for the spectral densities of the two-point functions for pseudoscalar and vector 
currents. The OPE in terms of the pole mass M /, calculated in [11] does not exhibit a perturbative hierarchy, therefore 
one rearrange the OPE in terms of the running MS-mass [12]. Then, the perturbative hierarchy of the correlation 
function starts to depend on tr; this feature allows one to choose the range of u where the perturbative hierarchy is 
visible. The negative effect of this rearrangement of the perturbative expansion is that, because of the truncation of the 
OPE series, the extracted decay constants acquire an unphysical dependence on the scale /i. In the charm sector this 
however does not lead to any serious problems. Figure 1 shows the dependence of the decay constants of the charmed 
pseudoscalar and vector mesons for the central values of all other OPE parameters after applying the algorithm for 
fixing the effective thresholds described above. One can see a weak /i-dependence of the decay constants of the 
pseudoscalar mesons mesons, whereas for vector mesons this /i-dependence is more pronounced. Averaging over the 
OPE parameters in their respective intervals and over the scale in the range 1 < ji [GeV] < 3 one arrives at the following 
results [19] 


fo = (208.3 ±7.3 0 PE±5 sy st) MeV, f Ds = (246.0± 15.7 OPE ± 5 syst ) MeV 

fir = (252.2±22.3 0 PE±4 sy st) MeV, f D * = (305.5 ±26.8 0 PE±5 S yst) MeV. (18) 

For the ratio we reported fo* / fo = 1.221 ± 0.080ope ± 0.008 sys t, which compares nicely with the lattice QCD result 
fo* / fD — 1.20 ±0.02. The results for the charmed mesons from other sum-rule analyses [17, 20] agree well with each 
other and with the results from lattice QCD [21], 
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FIGURE 1. Decay constants of D, D s D* and D* mesons depending on the scale ju. 
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FIGURE 2. Decay constants of B and B* mesons and the ratio f B * / f B depending on the scale fl. 


Beauty sector 

Similar to the charm sector, the OPE for pseudoscalar and vector currents containing the Z?-quark, does not show 
any perturbative hierarchy; there is no reason to assume that the unknown higher-order perturbative contributions are 
small. Rearranging the perturbative expansion in terms of the running mass introduces the dependence of the scale jj. 
and opens the possibility to choose the working range of fl in which the perturbative hierarchy is explicit thus allowing 
to hope the unknown higher orders do not contribute substantially to the correlation function. 

In the ^-sector one encounters two interesting features of the sum-rule analysis: 

• The sum-rule results for the beauty-meson decay constants correlate very strongly with the fc-quark mass [22] 

Sf B /f B ~ -%Sm b /tn h , (19) 

mf, = Til b (jn b ). Making use of m b = 4.18 ±0.03 GeV [23] leads to /« > 210 MeV, in clear tention with the recent lattice 
QCD results for f B ~ 190 MeV. Combining our sum-rule analysis with fn and f Bs from lattice QCD yields [22] 

m b = (4.247 ±0.027 ±0.01 l (syst) ) GeV. (20) 

The sum-rule results for the decay constants corresponding to this value of the fe-quark mass read 

f B = (192.0± 14.3 OP e± 3.0s y st) MeV, f B , = (228.0± 19.4 OPE ±4.0 sy st) MeV (21) 

• For the decay constant of B*, one observes an unexpectedly strong /i-dependence [24]: Averaging over the scale 
range 3 < ji [GeV] < 6 leads to 

/W/b = 0.923 ± 0.059, f B */f Bs = 0.932 ± 0.047. 

Taking into account only low-scale results for 2.5 < ju[GeV] < 3.5, yields f B */f B = 0.994±0.01. The sum-rule 
analysis [17] also gives indications that f B * /f B < 1 (ses Table II of [17]). Surprisingly, the QCD sum-rule prediction 
for f B * jfit is below the corresponding results from lattice QCD, which seem to favour a value slightly above unity 
[21, 25]. Clearly, such tension calls for further detailed investigations. 


/i -dependence of the physical quantities 

The heavy-light correltors are known with an impressive three-loop accuracy and are therefore rather weakly 
sensitive to the variations of the scale. Nevertheless, the dual correlator of the vector currents which includes the low- 
energy region of the Feynman diagrams only and, respectively, the vector-meson decay constants are rather sensitive 
to the choice of the scale. In many cases this scale-dependence is the main sourse of the OPE uncertainty in the decay 
constants. We should mention that in some publications the ^-dependence is treated in a specific way [20]: one just 
chooses one scale at which the decay constant has, e.g., an extremum in u, and provides the results for this very scale 
assigning no theoretical uncertainty to the scale fixing. This of course reduces strongly the total uncertainty of the 
decay constant obtained with the sum-rule technique but from our point of view such a treatment is not justified: the 
(unphysical) jj .-dependence is an effect of the truncation of the OPE series and thus reflects an essential feature of 
QCD. Any of the scale for which a reasonable perturbative hierarchy is seen, may be used for the determination of 
the hadron parameter; the unpleasant /i-dependence of the sum-rule results should be thus properly reflected in the 
theoretical uncertainty of the hadron parameter obtained using a QCD sum rule. 







SUM RULES FOR THREE-POINT VACUUM CORRELATION FUNCTIONS 


Let us now discuss the calculation of the meson elastic and transition form factors from the three-point vacuum 
correlation functions [26, 27], The basic object in this case has the form 

T{p 2 ,p' 2 ,q 2 ) = J (Q.\T(j(x)j(0)j(y))\Q.)exp(-ipx)exp(-ip'y)dxdy. (22) 

The three-point Green function in full QCD contains the double pole related to the mesons in the p 2 and //^-channels 
in the timelike region. The residue in this double pole is the form factor of interest. The Green function in the 
spacelike region may be calculated using the same method as the two-point function, i.e. by performing the OPE. 
One represents the Green function r(p 2 ,p' 2 ,q 2 ) as a double spectral integral in p 2 and p' 2 , performs the double Borel 
transform p 2 —> i and p' 2 —y x' (which, similar to the two-point function, kills the subtraction terms and suppresses 
the contributions of the excited states), equate to each other the OPE and the hadron representations for F(p 2 ,p' 2 ,q 2 ), 
and use duality property to isolate the ground-state contribution, thus relating the meson form factor to the low-energy 
region of the triangle diagrams of perturbative QCD and power corrections given through the condensates. For instance, 
the pion elastic form factor, in which case one sets x = X 1 , has the form [27] 


Seff(£? 2 ,T) we 2 ,r) 


i + s 2 T ( oc s G 2 ) 4na s (qq) 


i{Q 2 )fl= j j d^id52Apert(5i,52,6 2 )e ~^ T + X + 

0 0 

Apert {S1,S2,Q 2 ) = A (0) (.?[ , S 2 , Q 2 ) + a.vA (1) («t , S 2 , Q 2 ) + ' 
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(23) 


An essential feature of the three-point sum rule is that the effective threshold now depends on the Borel parameter x 
and the momentum transfer Q [28-30]; obviously, one faces a serious problem of finding appropriate algorithms for 
fixing ,v e ir(Q 2 . Tj. It should be understood that the effective continuum threshold for the form factor differs from the 
effective threshold for the decay constant. 1 

For large Q 2 , power corrections calculated in terms of the local condensates rise as polynomials with Q 2 , thus 
preventing a direct use of the sum rule (23) at large Q 2 . There are essentially only two possibilities to study the region 
of large Q 2 starting with the vacuum correlators; 

• use nonlocal condensates which are aimed at the resummation of the local condensate effects [31, 32], 

• work in the so-called local-duality (LD) limit x = 0 [31]. A specific feature of this limit is that all power corrections 
vanish in this limit and details of non-perturbative dynamics are hidden in one complicated object - the g 2 -dependent 
effective threshold s e ff(<2 2 )- 

A similar treatment may be performed for, e.g., the n° —> yy* transition form factor [33-35] for which one obtains 
the single spectral representation in the LD limit: 


%f(2 2 ) 

Fny(Q 2 )fn= j ds Op ert (s, Q 2 ) (24) 

o 

Due to properties of the spectral functions A pe rt(si A2>2 2 ) a ''d (Tpen A, Q 2 ), the form factors obey the factorization 
theorems 


Fn(Q 2 ) -> 8na s (Q 2 )f 2 /Q 2 , F ny (Q 2 ) -»• V2f n /Q 2 , f x = 130 MeV (25) 

as soon as the effective thresholds satisfy 

*eff(G 2 ^-)=^eff(e 2 ^-)=47r 2 / 2 . (26) 

Remarkably, due to the QCD factorization theorems for the hard form factors, the effective thresholds at Q 2 —> °° are 
given through the decay constants of the participating mesons. It should be emphasized that the only feature of theory 
relevant for this property of .v e ff(Q 2 ) is factorization of hard form factors. 


The effective thresholds for the baryon form factors are strongly sensitive to the choice of the interpolating current for a specific baryon. 



For finite Q 2 , the effective thresholds s e ff(0 2 ) and %f(2 2 ) depend on Q 2 and differ from each other [36, 37]. 
Nevertheless, setting s nne ff(Q 2 ) = s rme ff(Q 2 —> °°) for all not too small Q 2 [27] provides an approximate parameter- 
free prediction for the form factors which is becoming increasingly accurate as soon as Q 2 increases. The results of 
[37] give convincing evidences that s e ff(<2 2 ) and .s e ff (Q 2 ) are close to their asymptotic values already at relatively low 
values Q 2 ~ 4 — 8 GeV 2 . 

Thus, the LD approximation for the form factors—which requires as its crucial ingredient the knowledge of 0(1) 
and 0(a s ) double spectral densities—is increasingly accurate in the region not too close to zero recoil. The LD 
approximation is very promising for the application to, e.g., heavy-to-light weak form factors. A still missing ingredient 
here is the two-loop 0(a s ) double spectral density of the triangle diagram for different currents and arbitrary quark 
masses in the loop. This is a really challenging calculation which however opens the possibilities of very interesting 
applications. So far the only known results correspond to all massless quarks in the loop [38] and to HQET [39, 40], 


SUM RULES FOR THE EXOTIC POLYQUARK CURRENTS 

The OPE for the correlation functions of the exotic polyquark currents involving 4 (or 5) quark fields of the type 

D(x) = qi(x)Oq2(x)q3(x)Oq4(x) (27) 

where O is an appropriate combination of the Dirac matrices and possibly also of the (covariant) derivatives, have 
specific features compared to the OPE for the bilinear currents of the form j(x) = q \ (x) Oqo (x) used for usual 
“nonexotic” mesons. Namely, the lowest-order 0(1) contribution to the OPE for any correlator involving the exotic 
current, e.g. II dd = (0|r(D(x)D(0)|0), is given by the disconnected diagrams. As known from the general features 
of the Bethe-Salpeter equation and also emphasized recently by Weinberg [41], these disconnected diagrams are not 
related to the exotic bound states. The connected diagrams relevant for the exotic states emerge in the OPE for any 
correlator at the order O ( cx s ) and higher; therefore for the analysis of the exotic states the knowledge of the radiative 
corrections is mandatory. This makes the analysis of the exotic states a more technically involved problem than the 
analysis of the normal hadrons. 

Nevertheless, due to the fact that the observed exotic states are narrow, the procedure of extracting their parameters 
from the OPE has the same features and the same challenges as for the normal hadrons. Our experience in the analysis 
of the usual hadrons proves that a truncated OPE for the correlation function does not allow one to study at the same 
time both the existence of the isolated ground state and of its properties. However, if the mass of the narrow bound 
state is known, the method of sum rules allows one to obtain reliable predictions for its decay constants and the form 
factors. 


Structure of the exotic tetraquark states 

Obviously, the exotic tetraquark states may have a rather complicated “internal” structure; two most popular 
scenarios of this structure are a confined tetraquark state (i.e. a bound state in a confining potential between the 
two color-triplet diquarks) and a molecular “nuclear-physics like” bound state in the system of two colorless mesons. 

However, an important question about the structure of the exotic state—which to large extent determines also 
its production mechanism—is not easy to answer [42]: (i) by a combined color-spinor Fierz rearrangement of the 
tetraquark interpolating current D(x) one can write it either in diquark-antidiquark or meson-meson form; (ii) the 
same quantum numbers of the exotic interpolating current may be obtained by different combinations of its diquark- 
antidiquark or meson-meson bilinear parts. 

The simplest characteristic of a usual meson is its decay constant, i.e. the transition amplitude between the vacuum 
and the meson induced by its interpolating current; for a heavy quarkonium state the decay constant is analogous to its 
wave function at the origin y/(r = 0). 

For an exotic tetraquark state one should considers the connected self-energy functions 

n DD = (0|r(D(x)D(0)|0) = (DD) (28) 

and study the corresponding sum rules. However, for an exotic state one may obtain a set of the decay constants, 
related to different structure of the interpolating current with the quantum numbers of the exotic tetraquark of interest. 



The answer to the question of the dominant structure of the tetraquark may be given only by the analysis of a large set 
of the decay constants. 

• As the first step, one needs to study systematically the interpolating currents for tetraquark currents with different 
quantum nembers. As the next step one can calculate the set of YIdd- Because of the factorization property of the two- 
point function of the local tetraquark currents [43], the radiative corrections to II/ 5/5 are given via radiative corrections 
to the various two-point functions of the bilinear quark currents. For some of these two-point functions (namely, (VV) 
and (AA)) the radiative corrections are well-known, for some of them (such as (TT), T is the tensor bilinear current) 
these corrections should be calculated. 

• Then, the set of the sum rules for different two-point functions Udd should be studied and only then the answer 
about the structure of the observed narrow exotic candidates may be obtained. Especially interesting cases here are the 
narrow charged tetraquark Z (4430) (J r = 1 + and the width ~ 45 MeV, valence-quark content ccud) and X(3872) 
(J PC = 1 ++ X(3872), the width < 24 MeV). 

Another interesting possibility—so far not discussed in the literature—is considering nonlocal interpolating currents 
for the exotic mesons. The nonlocality of the interpolating currents should allow one to access in a better way subtle 
details of the tetraquark structure. 


Strong fall-apart decays of the exotic tetraquark states 


In the last decade, QCD sum rules have been extensively applied to the analysis of strong decays of exotic multiquark 
states (see e.g. [44, 45] and references therein). The basic object for the analysis of these decays in QCD is the three- 
point functions of the type 

T(p 1 p' 1 q) = J (0\T(D(0)j(xi)j(x 2 )\0)exp(-ip'xi-iqX2)dxidx2. (29) 

This correlator contains the triple-pole in the Minkowski region 
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(30) 


where dots stay for less singular terms. Here gxM t M-. is the three-hadron coupling which describes the X —> M\M 2 
transition; fx, fM t , and /m 2 are the decay constants of the mesons describing the strength of their their interaction 
with the interpolating current (X|D(0)|0) = fx and {M\ ,2171,2(0)|0) = f\ 2 (we omit here all Lorentz indices and for 
simplicity neglect the spins of the hadrons and of the interpolating currents). The OPE allows one to calculate the 
expansion of this correlator at the spacelike momenta far from the hadron thresholds. Again, the leading contribution 
in a s is given by a disconnected diagram (see Fig.3a) which factorizes and does not depend on the momentum of the 
exotic current p 2 at all: 


roPE (p 2 ,p' 2 ,< 7 2 ) = no' 2 )!!^ 2 ) + a s r connected (p 2 ,p' 2 ,q 2 ) (31) 




FIGURE 3. (a) The disconnected 0(1) diagram which does not depend on the variable p 2 relevant for the tetraquark properties; 

(b) One of the lowest-order connected 0(a s ) diagram which contributes to the tetraquark decay amplitude. 









Performing the Borel transform p 2 —> z, which comprises one of the steps of the sum-rule analysis, we see that the 
Borel image of the disconnected leading-order contribution vanishes. 2 Therefore any attempt to extract the tetraquark 
decay amplitude from the leading-order contribution is inconsistent. Relevant for the exotic-state properties are the 
0(a s ) corrections which are technically very difficult. This is a difficult calculation but it should be done before one 
may hope to get reliable predictions for the tetraquark properties. So far these corrections have been calculated only for 
the three-point function of the bilinear currents in two cases (i) for massless quarks and (ii) for infinitely heavy active 
quark and a massless spectator. For the 0(a s ) corrections to the three-point functions T, involving one tetraquark and 
two bilinear currents, no results exist in the literature. 

Nevertheless, the common feature of all previous calculations of these decays within QCD sum rules (e.g. [45, 46]) 
was the attempt to study the tetraquark (and pentaquark) decays basing on the factorizable leading-order contribution 
which intrinsically has no relationship with the tetraquark properties (which is clear both from the factorization prop¬ 
erty r{p,p',q) = n(// 2 )n(g 2 ) and from the large-A) behaviour of the QCD diagrams emphasized by Weinberg [41]. 
Therefore the existing analyses should be strongly revised by calculating and taking into account the nonfactorizable 
two-loop 0(a s ) corrections. 

Nonzero results based on the leading-order correlation function may be obtained only by a trick. Let us consider e.g. 
the decay Z —> t// + 7i~. One makes use of the tetraquark current j(x) = c{x)c{x)u{x)d(x) (we again omit the Dirac 
matrices for simplicity). The corresponding three-point correlation function of interest is 

r (p 2 ,p' 2 ,q 2 ) = J d 4 xd 4 yexp(—ip'x)exp(—iqy)(0\T(c(0)c(0)il(0)d(0),c(x)c(x),d(y)u(y))\0). (32) 

A nonzero result for the Borel transform of the disconnected zero-order contribution may be obtained by first 
considering the soft-pion limit q —7 0, i.e. p' = p. which gives for the disconnected contribution II(j? 2 )n(0) and 
then performing the Borel transform p 2 -X z. However, the decay rate obtained in this way is not really trustworthy. 

We therefore conclude that the “fall-apart” decay mechanism of exotic hadrons differs from the decay mechanism 
of the ordinary hadrons and requires the appropriate treatment within QCD sum rules. The calculation of the radiative 
corrections is mandatory for a reliable analysis of the properties of the exotic states. 


SUMMARY AND OUTLOOK 

In the recent years, great progress has been seen both in the calculations of the OPE series for various correlation 
functions and in the direction of formulating advanced algorithms for the extraction of the individual hadron parameters 
from these correlators. We could not discuss all the developments in this talk but let us try to mention in this summary 
the interesting open issues to be addressed in the future analyses: 

• Let us recall that combining moment QCD sum rules with experimental/lattice data gives the most accurate 
estimates of the heavy-quark masses [47]. 

• Hadron properties from 2-point functions: 

a. We have seen a visible progress in developing the new algorithms for extracting ground state parameters from 
the OPE of the correlators and gaining control over the systematic errors of the decay constants (finite-energy 
sum rules, Borel sum rules). Although it seems impossible to predict both masses and decay constants with a 
controlled accuracy, using the mass of the ground state as input, systematics can be controlled). 

b. We have encountered interesting puzzles in the / 7 -sector: 

(i) The Z>-quark mass 4.18 GeV [23] when used in the Borel sum rules for fg leads to tension with lattice results 
for f B - 

(ii) Unexpectedly strong scale-dependence of decay constants of vector mesons and of /«■ /fg even using the 
0(a 2 ) correlation function. 

c. Calculation of the decay constants of heavy-quarkonium states within the method of QCD sum rules is still not 
fully settled: The problem here is that the OPE for the doubly-heavy correlation functions contain relatively small 


2 We emphasize that for the decay of a usual hadron, the 0(1) contribution is given by a triangle diagram which depends on all three variables 
p 2 , p' 2 , f 2 and therefore of course does not vanish under the Borel transform p 2 —> T: for the decays of the usual hadrons the 0(1) contribution of 
the perturbative QCD indeed provides the dominant contribution to the decay of interest. 



nonperturbative power corrections. Therefore in QCD, the structure of OPE for the heavy-quarkonium system 
is somewhat similar to the structure of OPE for a purely coulomb system. Obviously, the algorithms adopted 
and tested for light or heavy-light hadrons in which cases the nonperturbative corrections are large, may work 
differently for heavy quarkonium states. This feature may be the origin of the tensions between the sum-rule 
predictions and the results from lattice QCD and other nonperturbative approaches for e.g. the decay constants of 
B c mesons and some charmonium states [20, 48]. A more critical analysis of the procedures of an isolation of the 
ground-state contribution from the correlation function and in particular of the way of obtaining the systematic 
uncertainties is necessary. 

d. Since the accuracy of the isolation of the ground-state contribution from the correlation function can be controlled, 
one may try to apply the method for the analysis of the excited states. Very little efforts in this direction have been 
done so far. 

• Meson elastic and transition form factors from three-point functions 

The Borel sum rules at x = 0 (the so-called local duality limit) open an interesting possibility of obtaining parameter- 
free predictions for the elastic and the transition form factors of light mesons in a broad range of the momentum 
transfers. The crucial ingredients necessary for these calculations are the 0(1) and 0(a s ) spectral densities of the 
triangle diagrams. As soon as these are known, the effective thresholds are determined in a unique way by the QCD 
factorization theorems for hard form factors. Assuming the effective thresholds to weakly depend on the momentum 
transfers—a hypothesis which finds support in the data for the pion form factors—one obtains the parameter-free 
predictions for the form factors in a broad range of the momentum transfers. Our analysis suggests that these 
representations for the form factors work with a few percent accuracy for Q 2 > a few GeV 2 . It seems very promising to 
apply the same ideas to heavy-to-light transition form factors. The main problem here is the necessity to calculate the 
radiative corrections to the triangle diagrams which is a very difficult task which needs serious efforts. As soon as this 
ambitious task is fulfilled, QCD sum rules could provide parameter-free predictions for the form factors, increasingly 
accurate with increase of Q 2 . 

• Baryon elastic and transition form factors 

The calculations for baryons are obviously technically extremely involved. Many sum-rule analyses of the baryon 
elastic and transition form factors have been presented in the recent years (see e.g. [49-52] and references therein). 
With a few exceptions (e.g. [49]), these calculations are based on the leading-order 0(1) correlation functions and use 
the traditional approaches to fix the effective thresholds, usually neglecting the T- and (^-dependence of the latter. 
These analyses are expected to provide reasonable ball-park estimates for the form factors; however, in most of the 
cases, the estimates of the OPE-errors (related to the uncertainties of the QCD parameters, to the missing radiative 
corrections, and in particular to a strong dependence on the scale p) and the systematic errors, related to the adopted 
procedures of fixing the effective continuum thresholds) are not done properly. Many efforts are still to be done in the 
domain of the baryon form factors. 

• Three-meson strong couplings of the type gD*Dx 

These quantities have been extensively addressed using three-point vacuum correlators and the corresponding sum 
rules. Again, the radiative corrections to the correlation functions have not been taken into account. Moreover, the 
results for the decay constants require extrapolations over large ranges of the momentum transfers. Therefore, one 
cannot expect good accuracy of these estimates. In many cases, the results from sum rules lead to an unrealistic 
picture of the SU (3)-breaking effects (see [53] and refs therein). For a real progress, one needs the calculation and the 
inclusion of the radiative corrections to the appropriate three-point functions. 

• Properties of the exotic tetraquark states 

The “fall-apart” decay mechanism of exotic hadrons differs from the decay mechanism of the ordinary hadrons and 
requires the appropriate treatment within QCD sum rules. In distinction to the decays of the usual hadrons, where the 
knowledge of the radiative corrections is necessary for improving the accuracy of the sum-rule form factor calculations, 
the calculation of the radiative corrections is mandatory for a reliable analysis of the properties of the exotic states. 
The decays of the exotic states are intrinsically unrelated to the 0(1) disconnected correlation functions; the results 
obtained from these 0(1) correlators cannot be treated as fundamental and reliable. 

From this summary of the recent advances and still open issues it seems obvious that the future progress in the 
sum-rule calculations of the properties of the usual and the exotic hadrons will be related (i) to the calculations of the 
radiative corrections to the correlation functions and (ii) to further development of the appropriate algorithms for the 
extraction of the properties of the individual hadrons from these correlators. 
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